Ion induced density bubble in a strongly correlated one dimensional gas 
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We consider a harmonically trapped Tonks-Girardeau gas of impenetrable bosons in the presence 
of a single embedded ion, which is assumed to be tightly confined in a RF trap. In an ultracold 
ion-atom collision the ion's charge induces an electric dipole moment in the atoms which leads to an 
attractive r~ 4 potential asymptotically. We treat the ion as a static deformation of the harmonic 
trap potential and model its short range interaction with the gas in the framework of quantum 
defect theory. The molecular bound states of the ionic potential are not populated due to the lack 
of any possible relaxation process in the Tonks-Girardeau regime. Armed with this knowledge we 
calculate the density profile of the gas in the presence of a central ionic impurity and show that a 
density bubble of the order of a micron occurs around the ion for typical experimental parameters. 
From these exact results we show that an ionic impurity in a Tonks gas can be described using a 
pseudopotential, allowing for significantly easier treatment. 
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Introduction - Since the days of superfluid helium 
the introduction of charged particles to a superfluid has 
proven to be a useful tool to explore fundamental physics. 
Statically they can lead to fundamental structures such 
as the celebrated ion snowballs and electron bubbles 
whereas studying the motion of an immersed ion can 
be used to characterise the superfluid nature of 4 He ■ 
Through the advances in the area of atom trapping and 
laser cooling the number of known superfluids and the 
environments in which they can be studied has in the 
last decade vastly increased. Combining these degener- 
ate fermionic and bosonic gases, which come in a wide 
range of geometries and dimensionalities Q with cooled 
and trapped single ions in Paul and Penning traps Q is 
therefore currently one of the most exciting and dynamic 
areas in physics. Such research is motivated by our need 
to understand the true quantum nature of matter, the 
quest for scalable quantum information processing and 
also the possibility of efficient quantum simulation Q. 
However, despite the achievements in both the areas of 
ion and atom physics, currently little is known about the 
physics of an ion interacting with a many-body atomic 
system in the ultracold regime @, 0| • 

By today most studies of atom-ion systems have fo- 
cused on the scattering properties of single atoms and 
single ions @-EIj- Controlled collisions of a single atom 
and a single ion have been studied in [llj with trap in- 
duced resonances predicted. Such resonances along with 
the strong atom ion interaction have shown to be useful 
for enhancing the speed and the fidelity of a collisional 
quantum gate [l2j ■ On the many-body level, it has been 
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proposed to use an ion as a scanning tunnelling micro- 
scope [lH which can in turn be used to measure the en- 
ergy distribution of, say, a Fermi gas in situ [lH • Despite 
these promising proposals there is still scarce knowledge 
of what is to be expected when an ion is introduced to an 
ultracold many-body system. As the physics of a many- 
body system is highly non-trivial even without the pres- 
ence of an ion, such studies are of paramount importance 
to shed light on the first results of the ongoing exper- 
imental progress [15l4l7l |. Two situations have recently 
been explored. In [6| the authors consider an ion in a ho- 
mogeneous condensate and study the capture of atoms 
into the bound states of the molecular potential. They 
predict the possibility of observing mesoscopic molecular 
ions, where hundreds of the atoms can become bound in 
molecular states of the ion-atom potential. On the other 
hand, in Q the authors consider a situation where the 
bound states are not accessible. They then estimate by 
using a thermodynamical as well as a microscopic argu- 
ment the number of condensate atoms associated with an 
ion. Both studies suggest strong analogies with the snow- 
ball states observed in the case of superfluid helium. In 
contrast to both of these studies, in this work we demon- 
strate that in the case of a harmonically trapped and 
strongly correlated one-dimensional gas, a micron sized 
bubble in the density may be observed around the ion. 

Strongly correlated one-dimensional quantum gases 
have recently become available in ultracold physics labs 
worldwide. One of the stand out achievements for cold 
atoms in low dimensions has been the first observation of 
the Tonks-Girardeau (TG) regime [Hdl]. In this limit 
the gas is defined to be a one-dimensional, strongly cor- 
related gas [2(|. As it is exactly solvable through the 
Bose-Fermi mapping theorem [20j it is very well suited 
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FIG. 1: Schematic showing (a) an positively charged ion em- 
bedded in a harmonically trapped Tonks-Girardeau gas and 
(b) the interaction potential for the ion-atom interaction un- 
der the Born-Oppenheimer approximation (dashed). The po- 
tential minimum is at r m in and the polarisation length of the 
ion is indicated by R* . The potential curve goes as r~ 4 at 
large distances. 



to study its behaviour in the presence of impurities. The 
physics of impurities interacting with a many-body sys- 
tem is of great importance in understanding many com- 
plex condensed matter phenomena and as ultra-cold sys- 
tems allow to tailor the interaction they have proven to 
be versatile model s yste ms for a large number of differ- 
ent areas in physics |2l| . In fact, a recent experiment in 
Cambridge has observed the motion of a spin impurity 
in a TG gas and reported interesting and complex dy- 
namics [22[. In this letter we add to the existing studies 
@j of ionic impurities interacting with atomic many- 
body systems by considering an ion embedded in a TG 
gas (see Fig. (Ha)). Utilising the Fermi-Bose mapping 
theorem [2(| we calculate the single particle density of a 
TG gas in the presence of an ion. We treat the ion's pres- 
ence as a perturbation on the harmonic trap and model 
its short range interaction in the framework of quantum 
defect theory. 

Ultracold atom-ion interactions - Under the Born- 
Oppcnhcimer approximation, the interaction between an 
atom and a positively charged ion, at large distance r and 
low energies, is given by the potential V(r) = — ae 2 /2r 4 . 
Here a is the dipolar polarizability of the atomic species 
and e is the charge on the electron. This approximation 
is valid only if the atom is in an atomic s state as for 
higher electronic states the Born-Oppenheimer potential 
significantly changes shape. A characteristic lengthscale 
of this interaction can be found by equating the interac- 
tion energy and the relative kinetic energy of an atom- 
ion pair. This gives the so-called polarisation length 



R* = y/2fiae 2 /ft 2 , where fi is the reduced mass and 
which allows to define the interaction energy of a typ- 
ical atom ion collision as E* = ft 2 /2/i(R*) 2 . 

The Hamiltonian for the axial dynamics of an atom 
and an ion in seperate one dimensional harmonic traps 
maybe written as (Tl| . 
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where V, n t is the full one dimensional interaction poten- 
tial, which at large distances has the same power law as 
in three dimensions i.e. V(|a;|) = —ae 2 /2x 4 . In many 
realistic experimental scenarios the frequency of the ion 
trap will be many orders of magnitude greater than that 
of the atom trap and it is reasonable to approximate the 
ion to be located quasi-statically at the centre of the atom 
trap. In this approximation one may neglect the kinetic 
and potential energy of the ion in (Q]). This allows us 
to write a single particle Hamiltonian which describes an 
atom in a harmonic trap which is deformed at the trap 
centre by the ionic interaction potential, 
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The solution of the Schrodinger equation for the single 
particle states is complicated by the fact that the real 
interaction potential Vint begins to deviate from the long 
range asymptotic law and begins to quickly diverge to- 
wards +00 (see Fig. QXb)). However, at this distance the 
effect of the harmonic trapping potential can safely be 
neglected and at short range we may approximate the 
scattering behaviour to be equivalent to a free atom-ion 
scattering event. The Schrodinger equation, in this re- 
gion, maybe written as 



where \x is the reduced mass and <p e and <p are the un- 
known scattering phases to be determined (see eqs. ^ 
and ([3]) below ). In one dimension the solutions for eq. ([3]) 
have even and odd parity. The asymptotic solutions, 
\x\ — > 0, to cq. ([3]) are given by fill] 
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The principle idea of the quantum defect theory is to 
replace the real interaction potential in eq. ©, Vint, 
with the asymptotic potential and use eqs. (01 and ([5]) 
as boundary condition for the numerical solution of the 
Schrodinger equation for Hamiltonian cq. © . In analogy 
to the three-dimensional case the short range phase may 
be related to the one-dimensional even (s— wave) and odd 
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FIG. 2: The first two trap states of a harmonic trap with a 
centrally embedded ion for a 135 Ba + - 87 Rb system for axial 
trapping frequency, ui — 70Hz. The inset for the n — 
state shows a zoom in of the region around zero where the 
wavefunction is rapidly oscillating. The short range phases 



are chosen to be 



and 



(p— wave) scattering lengths of an ultracold atom-ion col- 
lision via afn = li m fc->o[tan(0 eiO (fc))/fc] = — cot(0 e]O ). 
Unfortunately the exact values for the scattering lengths 
of current experimental systems are not yet known and 
we therefore take </> ej0 to be adjustable parameters for 
now. The set of single particle states obtained using this 
approach are the ones which are used to calculate the 
groundstate of a Tonks-Girardeau gas in the presence of 
a centrally embedded ion. 

gas - The effective one- 
of a gas of N bosons in 
can be written as Ti = 



The Tonks-Girardeau 
dimensional Hamiltonian 
trapping potential V(x) 
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12i<jV(\xi-Xj\), where m 

is the mass of a single atom. For low densities two body 
collisions dominate the dynamics, and for low tempera- 
ture we can approximate these by employing a point-like 
potential of the form V(|a;j — Xj |) = giD$(\xi — Here 
g±D is the one dimensional coupling constant, which is 
related to the three dimensional scattering length, a3D, 
via gin = 4 ^°^ p (l± — Ca3£)) _1 , where C is a constant 
of value C = 1.4603 . . . [23[ and l± is the transverse 
trapping length. In the Tonks-Girardeau limit of in- 
finitely strong repulsion (giD ~> oo), the problem may 
be mapped to the problem of non-interacting spinless 
fermions, which is the so called Fermi Bose mapping the- 
orem [20( ■ The main idea of the mapping is that one can 
treat the atom-atom interaction term by replacing it with 
the following boundary condition on the allowed bosonic 
wave- function, ^Sb{xi . . . xn) = if \xt — Xj\ =0 for 
i =/= j and 1 < i < j < N. As this is formally equivalent 
to the Pauli exclusion principle, one can solve for the as- 
sociated ideal fermionic wave- function 1 F f{xi, ■ ■ ■ , %n) = 
^= det ni j—i[ip n (xj)], and calculate the bosonic solu- 
tion from this by appropriate symmetrization, ^b = 
A(xi, . . . , xn)^f(xi, X2, ■ ■ ■ , Xn), where the unit anti- 
symmetric function is given by A = rii<i<j<jv s S a ( x i ~ 
Xj) [2fj| . As a consequence of this mapping the single 
particle density of a TG gas has the particularly simple 
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and thus all one needs to know to solve the many particle 
problem are the eigenstates of the single particle problem. 

Ion in a Tonks- Girardeau gas - By combining the ideas 
of the previous two sections we are now in a position to 
calculate the effect of a centrally embedded ion on the 
ground state of a harmonically trapped TG gas. The 
Hamiltonian of the combined system can be written as 

N / d 2 1 
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We now focus on obtaining the eigenstates of the sin- 
gle particle part of the Hamiltonian, i.e. eq. ((2]). Using 
the prescription of the Fermi-Bose mapping we neglect 
the interaction part of the Hamiltonian and numerically 
solve the single particle Schrodinger equation to find the 
eigenstates ipn and eigenvalues e n . For this we use the 
iterative Numerov method [HI and impose the eqs. (|4]) 
and |[SJ) as short range boundary conditions. 

Due to the infinitely diverging nature of the ionic po- 
tential, there exists a spectrum of bound molecular states 
between the ion and the atoms Q . In order for these to be 
populated, however, a relaxation process must be present 
that allows to take away the excess energy of the bind- 
ing process. In an ultracold gases this process would be 
equivalent to a three body recombination. However for 
the TG gas the density-density correlation function be- 
comes suppressed on the length scale of the interparticle 
distance [25j - and thus no relaxation process for capture 
is available. We therefore do not need to consider the 
bound state spectrum in the following. 

Our calculations are carried out for the experimentally 
relevant systems of 135 Ba+- 87 Rb and 174 Yb+- 87 Rb. We 
assume a typical axial trap frequency of to = 70Hz for 
a one dimensional gas of 87 Rb. The 135 Ba + and 174 Yb + 
ions have typical polarisation lengths of R* = 5544ao and 
6294ao, respectively, where ao is the Bohr radius. For the 
135 Ba + - 87 Rb system we show the first two single particle 
eigenstates in Fig. @ . One can see that the centre of the 
eigenstates is perturbed by the presence of the ion and a 
close up of this central region is shown in the inset for the 
ii = state. One can also see that in this region the wave 
function is rapidly oscillating, with the exact details of 
the oscillations depending on the short range phase. In 
the limit of such large trap lengths the effects of the oscil- 
lations on the physics are only marginal and only in the 
opposite limit to ours, if the traps were tight, do the os- 
cillations and hence the exact short range phase become 
important. The fact that the anti-symmetric states are 
affected means that a certain amount of p-wave scatter- 
ing is to be expected, however it will be very small in 
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FIG. 3: The single particle densities of a TG gas of twenty 
particles in the presence of a central ion for 135 Ba + - S7 Rb and 
174 Yb + - S7 Rb systems with a typical trapping frequency, u — 
70Hz (thick black line). The short range phases in each case 
are chosen to be <f> e = j and 4> = —j. The thin red line in the 
plots represents the result of a pseudopotential approximation 
for the ion, using a large value for the scattering length. 

our wide trap. With the single particle wavefunctions at 
hand, we can now calculate the density distribution of 
a many-particle gas using eq. (JB]). In Fig. ([3]) we show 
the density for the two separate systems of 135 Ba + - 87 Rb 
and 174 Yb + - 87 Rb. The thick black line shows the re- 
sults using the quantum defect theory and the Numerov 
method including the full atom-ion polarisation poten- 
tial. It can be seen that the ion's presence drastically 
perturbs the density distribution of the gas by causing 
a bubble in the density around its position at x = 0. 
The size of this bubble is of the order of lfim. At first 
this may seem counter-intuitive for an attractive impu- 
rity but one should recall that scattering is repulsive at 
short distances (see Fig. [T}d) and with the lack of access 
to the bound states of the molecular potential in the TG 
regime - we see no build up of atomic density around the 
ion. 

Let us finally note the possibility of treating atom-ion 
systems in the Tonks limit with a simpler approach, but 
approximating the exact interaction with a point-like po- 
tential [2(| . The polarisation lengths for both 135 Ba + and 



174 Yb + are, 0.29/ito and 0.33/im, and therefore smaller 
than trap lengths for one dimensional systems, which for 
the traps considered here are ~ 1.3/im. This suggests 
that a localised pseudopotential approximation can be 
used and in Fig. [3]we show this. In Fig. [3] the pseudopo- 
tential approximation is shown by the thin red line and 
one can see that the two descriptions are in good agree- 
ment despite the fact that the scattering length in this 
case is not much smaller than the trap length. There- 
fore one can expect that the pseudopotential approxi- 
mation will encapsulate much of the essential physics on 
the many-body level. In particular, the analytic accessi- 
bility of the pseudopotential model allows one to easily 
calculate the correlation functions and momentum dis- 
tribution 12711 of a TG gas in the presence of a localised 
impurity |26| . The pseudopotential model accurately de- 
scribes the localised nature of the disturbance but due 
to it's zero width nature it underestimates the size of 
the bubble. It is also important to mention, where den- 
sities are concerned, that these results are also true for 
spin polarised samples of non-interacting fermions in one 
dimension. 

Conclusions - In this work we have shown that the 
presence of a single ion in a TG gas leads to the formation 
of a density bubble around the impurity. Using quantum 
defect theory we have calculated the exact eigenstates of 
an atom-ion system in one dimension and then applied 
the Bosc-Fcrmi mapping theorem to calculate the many- 
particle density. The bubble was found to be of the order 
of a /im, and should therefore be observable in typical 
experiments. We also outlined that the density can be 
very well approximated using a pseudopotential model 
under typical experimental constraints, which will allow 
to connect the atom-ion systems in the Tonks limit to a 
large number of results already existing in the literature. 
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